Neutrino-induced nucleosynthesis of A > 64 nuclei: The z/p-process 
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We present a new nucleosynthesis process, that we denote i^p-process, which occurs in super- 
novae (and possibly gamma-ray bursts) when strong neutrino fluxes create proton-rich ejecta. In 
this process, antineutrino absorptions in the proton-rich environment produce neutrons that are 
immediately captured by neutron-deficient nuclei. This allows for the nucleosynthesis of nuclei with 
mass numbers A > 64. Making this process a possible candidate to explain the origin of the solar 
abundances of ^^■^''Mo and ^'^'^^^Ru. This process also offers a natural explanation for the large 
abundance of Sr seen in an hyper-metal-poor star. 

PACS numbers: 26.30.+k, 97.60.Bw 



Supernova explosions (the cataclysmic endpoint of stel- 
lar evolution) produce iron and neighboring nuclei, un- 
derlined by their lightcurves powered by radioactive de- 
cay of ^^Ni. The production of elements beyond Fe has 
long been postulated by three classical processes, the r- 
and the s-process (caused by rapid or slow neutron cap- 
ture) and the p-process, standing either for proton cap- 
ture or alternative means to produce heavy neutron de- 
ficient, stable isotopes [HQl- The s-process acts during 
stellar evolution via neutron captures on Fe produced 
in previous stellar generations (thus being a "secondary 
process"). The location and/or operation and uniqueness 
of the r- and p-process in astrophysical sites are still a 
subject of debate. The r-process is required to be a pri- 
mary process in stellar explosions Q, meaning that the 
production of such elements is independent of the ini- 
tial heavy element content in the star. Recent galactic 
chemical evolution studies of Sr, Y, and Zr 4j suggest 
the existence of a primary process, denoted "lighter ele- 
ment primaryprocess" (LEPP), that is independent of 
the r-process |^ and operates very early in the Galaxy. 
Most of the p-nuclei are thought to be produced in hot 
(supernova) environments, where disintegration of pre- 
existing heavy elements (thus being also a secondary pro- 
cess) due to black-body radiation photons can account 
for the heavy p-nuclei but underproduces the light ones 
(see e.g. ref. [a Bi Currently, the mechanism for 

the production of the light p-nuclei, ^^'^''Mo and ^^'^'^Ru, 
is unknown, however chemical evolution studies of the 
cosmochronometer nucleus ^^Nb , imply a primary su- 
pernova origin for these light p-nuclei. 

Observations of extremely "metal-poor" stars in the 
Milky Way provide us with information about the nucle- 
osynthesis processes operating at the earliest times in the 



evolution of our Galaxy. They are thus probing super- 
nova events from the earliest massive stars, the fastest 
evolving stellar species. The recently discovered hyper- 
metal-poor stars in the Milky Way [lOl llll | may witness 
chemical enrichment by the first generation of faint mas- 
sive supernovae which experience extensive matter mix- 
ing (due to instabilities) and fallback of matter after the 
explosion However, the detection of Sr/Fe, exceed- 
ing 10 times the solar ratio, in the most metal-poor star 
known to date I l0| suggests the existence of a primary 
process, producing elements beyond Fe and Zn. 

In this Letter, we present a new nucleosynthesis pro- 
cess that will occur in all core-collapse supernovae and 
could explain the existence of Sr and other elements be- 
yond Fe in the very early stage of galactic evolution. We 
denote this process "z/p-process" and suggest it as a can- 
didate for the postulated lighter element primary process 
LEPP 01 . It can also contribute to the nucleosynthesis of 
light p-process nuclei. Here, we consider only the inner 
ejecta of core-collapse supernovae, but the winds from 
the accretion disk in the coUapsar model of gamma-ray 
bursts 0, Q, 0, 01 may also be a relevant site for 
the i/p-process. The z^p-process is primary and is asso- 
ciated with explosive scenarios. It occurs when strong 
neutrino fluxes create proton-rich ejecta. After the pro- 
duction of Fe-group elements, continued antineutrino ab- 
sorptions by free protons produce free neutrons subject 
to immediate capture on neutron-deficient nuclei with 
small proton-capture cross sections and long beta-decay 
half-lives. This is distinct from earlier suggestions for the 
production of p-nuclei in the so-called neutrino wind that 
develops in later phases of a supernova explosion [l7lll8| . 
because the captured neutrons are directly delivered by 
antineutrino absorption on free protons in a proton-rich 
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environment. 

As a full understanding of the core collapse super- 
nova mechanism is still pending and successful explo- 
sion simulations are difficult to obtain the compo- 
sition of the innermost ejecta - directly linked to the 
explosion mechanism - remained to a large extent un- 
explored. Recent supernova simulations with accurate 
neutrino transport [i^, 0, 112 show the presence of 
proton-rich neutrino-heated matter, both in the inner 
ejecta |^ and in the early neutrino wind from the 
proto- neutron star j^^. This matter is subject to a large 
neutrino energy deposition by the absorption of neutrinos 
and antineutrinos with initially similar intensities and en- 
ergy spectra. As soon as the heating lifts the electron de- 
generacy, the reactions Ve+n «^ p+e~ and Ve+P ^ ?i+e+ 
drive the composition proton-rich due to the smaller mass 
of the proton |23ll2^ . (n, p, e", e"*", Ve, f'e denote the neu- 
tron, proton, electron, positron, neutrino, and antineu- 
trino respectively). Such proton-rich matter with Y^, the 
number of electrons or protons per nucleon, larger than 
0.5 will always be present in core-collapse supernovae 
explosions with ejected matter irradiated by a strong 
neutrino flux, independently of the details of the explo- 
sion [i^l . As this proton-rich matter expands and cools, 
nuclei can form resulting in a composition dominated 
hy N = Z nuclei, mainly ^^Ni and ""^He, and protons. 
Without the further inclusion of neutrino and antineu- 
trino reactions the composition of this matter will finally 
consist of protons, alpha-particles, and heavy (Fe-group) 
nuclei (in nucleosynthesis terms a proton- and alpha-rich 
freeze-out), with enhanced abundances of ''^Sc, ''^Ti, and 
®^Zn [ISIijI- these calculations the matter flow stops 
at ^"^Ge with a small proton capture probability and a 
beta-decay half-life (64 s) that is much longer than the 
expansion time scale (~ 10 s) |2^ . 

Synthesis of nuclei heavier than ^ = 64 is possi- 
ble in proton rich ejecta if the entropy per nucleon is 
in the range s « 150-170 ks (where fc^ is the Boltz- 
mann constant) Such large entropies are, how- 

ever, not attained in core-collapse supernovae simulations 
with detailed neutrino transport which give s « 50- 
75 ks HE Here we show that the synthesis of 

nuclei with A > 64 can also be obtained with realistic 
entropies, if one explores the previously neglected effect 
of neutrino interactions on the nucleosynthesis of heavy 
nuclei. When interactions with neutrinos and antineu- 
trinos are considered for both free and bound nucleons 
the situation becomes dramatically different. N ^ Z 
nuclei are practically inert to neutrino capture (convert- 
ing a neutron in a proton), because such reactions are 
endoergic for neutron-deficient nuclei located away from 
the valley of stability. The situation is different for an- 
tineutrinos that are captured in a typical time of a few 
seconds, both on protons and nuclei, at the distances 
at which nuclei form (~ 1000 km). This time is much 
shorter than the beta-decay half-life of the most abun- 
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FIG. 1: Evolution of the abundance of neutrons, protons, 
alpha-particles, and ^^Ni in a nucleosynthesis trajectory re- 
sulting from model B07 of reference |23|. The abundance, Y , 
is defined as the number of nuclei of the species i present di- 
vided by the total number of nucleons which is conserved dur- 
ing the calculation. The solid (dashed) lines display the nucle- 
osynthesis results which include (omit) neutrino and antineu- 
trino absorption interactions after nuclei are formed. The 
abscissa measures the time since the onset of the supernova 
explosion. 

dant heavy nuclei reached without neutrino interactions 
(e.g. ^^Ni, 64Ge). As protons are more abundant than 
heavy nuclei, antineutrino capture occurs predominantly 
on protons, causing a residual density of free neutrons of 
lO^'^-lO^^ cm~'^ for several seconds, when the tempera- 
tures are in the range 1-3 GK. This effect is clearly seen 
in figure ^ where the time evolution of the abundances 
of protons, neutrons, alpha-particles and ^^Ni is shown 
(^^Ni serves to illustrate when nuclei are formed). The 
dashed lines shows the results for a calculation where 
neutrino absorptions are neglected once the temperature 
drops below 6 GK. This allows to study the effect of 
neutrino absorptions in the latter phase of nucleosynthe- 
sis when the i^p-process acts without changing the initial 
phases where Ye is determined. Without the inclusion 
of antineutrino capture the neutron abundance soon be- 
comes too small to allow for any capture on heavy nuclei. 
The figure also compares the evolution of Y^ . 

In our studies we use the detailed neutrino spectral in- 
formation provided by neutrino radiation hydrodynam- 
ical calculations to determine the neutrino antineutrino 
absoption rates at each point of the nucleosynthesis tra- 
jectory (temperature, density and radius). Our network 
calculations follow the detailed abundances of 1435 iso- 
topes between Z = 1 and Z = 54, which allows an ac- 
curate treatment of the changes in composition induced 
by neutrino interactions. However, our network calcu- 
lations follow the Ye evolution of the hydrodynamical 
calculations only till the moment when alpha particles 
form. At this time, the determination of the Yf. value in 
the hydrodynamical studies is plagued by an error in the 
Lattimer-Swesty equation of state |26ll27l|. which we had 
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FIG. 2: The panel shows the isotopic abundances for model 
B07 of reference relative to solar abundances |29|l . The 
filled circles represent calculations where (anti)neutrino ab- 
sorption reactions are included in the nucleosynthesis while 
for the open circles neutrino interactions are neglected. 



adopted in the hydrodynamical calculations. This error 
results in an underproduction of alpha particles which 
suppresses the occurence of an alpha-effect which drives 
the Ye value closer to 0.5 While such alpha ef- 

fect does not occur in the hydrodynamical calculations as 
the computed alpha abundance is too low, it is present 
in the network calculations. However, we stress that, in 
contrast to expectation 0| , the alpha effect is no obsta- 
cle for the nucleosynthesis of heavy nuclei, because once 
heavy nuclei form the neutrons are captured by heavy 
neutron-deficient nuclei instead of forming deuterium and 
later alpha particles. 

The neutrons produced via antineutrino absorption 
on protons can easily be captured by neutron-deficient 
N Z nuclei (for example ^^Ge), which have large neu- 
tron capture cross sections. The amount of nuclei with 
A > produced is then directly proportional to the 
number of antineutrinos captured. While proton cap- 
ture, (p, 7), on ^''Ge takes too long, the {n,p) reaction 
dominates (with a lifetime of 0.25 s at a temperature of 
2 GK), permitting the matter flow to continue to nuclei 
heavier than ^''Ge via subsequent proton captures with 
freeze-out at temperatures around I GK. This is differ- 
ent to r-process environments with < 0.5, i.e. neutron- 
rich ejecta, where neutrino capture on neutrons provides 
protons that interact mainly with the existing neutrons, 
producing alpha-particles and light nuclei. Their cap- 
ture by heavy nuclei is suppressed because of the large 
Coulomb barriers 0, H^. Consequently, in r-process 
studies an enhanced formation of the heaviest nuclei does 
not take place when neutrino interactions are included. 
In proton-rich ejecta, however, antineutrino absorption 
produces neutrons that do not suffer from Coulomb bar- 
riers and are captured preferentially by heavy neutron- 
deficient nuclei. 



Figure|21shows the composition of supernova ejecta ob- 
tained with the hydrodynamical model B07 described in 
detail in ref. [l^. In addition to the proton-rich condi- 
tions in the innermost ejected zones visible in simulations 
by different groups in EH 113, our models consistently 
include neutrino-absorption reactions in the nucleosyn- 
thesis calculations allowing for the occurrence of the vp- 
process. However, in our stratified spherically symmetric 
models the accretion rate is rapidly reduced (and with 
this the neutrino luminosities) with the onset of the ex- 
plosion. In a more realistic scenario considering convec- 
tive turnover in the hot mantle, continued accretion is 
expected to maintain a large neutrino luminosity beyond 
the onset of the explosion and to further support the i^p- 
process 

In order to understand the sensitivity of our results 
one must consider the dependence of the I'p-process on 
the conditions during the ejection of matter in super- 
nova explosions. There are several essential parameters 
in addition to the entropy s. One is the Kg-value of the 
matter when nuclei are formed. The larger the Kg-value, 
the larger is the proton abundance, producing a larger 
neutron abundance for the same antineutrino flux during 
the i^p-process. This permits a more efficient bridging 
of beta-decay waiting points by (n,p)-reactions in the 
flow of proton captures to heavier nuclei. The location 
(radius r) of matter during the formation of nuclei and 
the ejection velocity also influence the z^p-process by de- 
termining the intensity and duration of the antineutrino 
flux which the matter will experience. A location closer 
to the surface of the protoneutron star and/or a slow 
ejection velocity leads to an extended antineutrino ex- 
posure that allows for an increased production of heavy 
elements. Finally, the long-term evolution of the neutrino 
luminosities and energy spectra during the cooling phase 
of the proto-neutron star plays an important role. These 
factors are poorly known due to existing uncertainties in 
the supernova explosion mechanism. 

To test the dependence of the nucleosynthesis on these 
parameters we have also carried out parametric calcula- 
tions based on adiabatic expansions similar to those used 
in refs. HO] but for a constant realistic entropy per 
nucleon s = 50 fcs. This allows exploration of the sensi- 
tivity of the nucleosynthesis without the need to perform 
full radiation-hydrodynamical calculations. An example 
is given in figure |31 which shows the dependence of the 
p-nuclei abundances as a function of the Y^, value of the 
ejected matter. The different Yg values have been ob- 
tained by varying the temperatures of the neutrino and 
antineutrino spectra assuming Fermi-Dirac distributions 
for both. Close to Y^ = 0.5 (and below) essentially no 
nuclei beyond A — G-i are produced. Nuclei heavier than 
A — 64: are only produced for Y,, > 0.5, showing a very 
strong dependence on Yg in the range 0.5-0.6. A clear 
increase in the production of the light p-nuclei, ^^^^^Mo 
and ^^^^*Ru, is observed as Ye gets larger. However, addi- 
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FIG. 3: Light p- nuclei abundances in comparison to solar 
abundances as a function of Y^. The Ve- values given are the 
ones obtained at a temperature of 3 GK that corresponds to 
the moment when nuclei are just formed and the i^p-process 
starts to act. 

tional work is necessary, including an exploration of the 
nuclear physics uncertainties, to decide whether or not 
the solar system abundances of light p-nuclei are due to 
the t'p-process. We obtain similar variations (not shown 
in the figure) in the production of Sr, Y and Zr, in agree- 
ment with the large scatter of the abundances seen in low 
metallicity stars 

After we have presented our results in several confer- 
ences, the j/p-process has also been explored by Pruet 
et al |3l| . who confirm our findings. They have con- 
sidered the nucleosynthesis in the early proton-rich neu- 
trino wind in a supernova explosion. The conditions- 
there, larger entropies {'^ 75 kb) and nuclei forming at 
shorter distances (several hundreds km) from the neutron 
star, extend the process described here to even heavier- 
nuclei. 

All core-collapse supernova explosions, independent of 
existing model uncertainties, will eject hot, explosively 
processed matter subject to strong neutrino irradiation. 
We argue that in all cases the ^-p-process will operate 
in the innermost ejected layers producing neutron defi- 
cient nuclei above A = 64. As the innermost ejecta, this 
matter is most sensitive to the details of individual ex- 
plosions, thus their abundances will vary noticeable from 
supernova to supernova (e.g. as a function of stellar mass, 
rotation, etc.). The final amount of matter ejected will 
also depend on the intensity of the fallback, but as dis- 
cussed in ref . ^2 > mixing before fallback will always lead 
to the ejection of elements synthesized, even in the in- 
nermost layers. Ref. explains the abundances seen 
in HMP stars by the ejecta of faint/weak core-collapse su- 
pernovae. Such faint supernova will generally have small 
expansion velocities favoring an enhanced production of 
j/p-elements, offering an explanation for the presence of 
Sr in the star HE 1327-2326 Our studies of the 



fp-process show that the elements between Zn and Sr 
should be co-produced together with Sr. The observa- 
tion of these elements, which with the exception of Ge 
and Rb are not detectable from the ground in optical 
lines, but possible from space in the infrared or near ul- 
traviolet (e.g. the Hubble Space and Spitzer Space tele- 
scopes), can provide support for the occurrence of the 
i/p-process at early times in the Galaxy. Such observa- 
tions can contribute valuable information about the con- 
ditions experienced by the inner supernova ejecta in order 
to constrain current theoretical models of supernova ex- 
plosions. Further studies are required to fully understand 
the ^-p-process contribution to the chemical evolution of 
the Galaxy. In this paper, we have discussed the inner- 
most proton-rich supernova ejecta before the emergence 
of the neutrino-wind from the protoncutron star. This 
neutrino-wind will initially be proton-rich [24l | but will 
turn neutron-rich in its later phases allowing for the syn- 
thesis of r-process nuclei 0, 0, The variations 
in the contribution of the t'p-process (represented by Sr, 
Y, and Zr) and the r-process (producing the heaviest el- 
ements up to Th and U) 0, can shed light on the 
connection of both of these processes and provide infor- 
mation about the class of supernovae that produce the 
heavy r-process nuclei. 
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